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r ~7 l ! Abstract 

In this paper we will consider Rohde's general form of {l}-inverse of a matrix 

A. The necessary and sufficient condition for consistency of a linear system 

Ax = c will be represented. We will also be concerned with the minimal 

number of free parameters in Penrose's formula x = A^c + (I — A^A)y 

for obtaining the general solution of the linear system. This results will 

be applied for finding the general solution of various homogenous and non- 

-homogenous linear systems as well as for different types of matrix equations. 
O 
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CO ' 1. Introduction 



In this paper we consider non-homogeneous linear system in n variables 

Ax = c, (1) 

where A is an m x n matrix over the field C of rank a and c is an m x 1 
matrix over C. The set of all m x n matrices over the complex field C will be 
denoted by C mxn , m,nGN. The set of all m x n matrices over the complex 
field C of rank a will be denoted by C™ xn . For simplicity of notation, we will 
write Aj_> (Alj) f° r the i th row (the j column) of the matrix A G C mxn . 
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Any matrix X satisfying the equality AX A = A is called {l}-inverse of A 
and is denoted by A^\ The set of all {l}-inverses of the matrix A is denoted 
by A{1}. It can be shown that ^4{1} is not empty. If the n x n matrix A 
is invertible, then the equation AX A = A has exactly one solution A -1 , so 
the only {l}-inverse of the matrix A is its inverse A -1 , i.e. v4{l}={y4 -1 }. 
Otherwise, {l}-inverse of the matrix A is not uniquely determined. For 
more informations about {l}-inverses and various generalized inverses we 
recommend A. Ben-Israel and T.N.E. Greville [l| and S.L. Campbell and CD. 
Meyer Q. 

For each matrix A e c™ xn there are regular matrices P G C nxn and 
Q E C mxm such that 



QAP = E a 



\h 


" 









(2) 



where I a is a x a identity matrix. It can be easily seen that every {l}-inverse 
of the matrix A can be represented in the form 



A® = P 



\I« 


u ' 


V 


w 



Q 



(3) 



where U = [uij], V = [vy] and W = [ivy] are arbitrary matrices of corre- 
sponding dimensions a x (m — a), (n — a) x a and (n — a) x_(m — a) with 
mutually independent entries, see C. Rohde [8| and V. Peric J7|. 

We will generalize the results of N. S. Urquhart [9J. Firstly, we explore 
the minimal numbers of free parameters in Penrose's formula 



x = A {1) c+(I-A {1) A) 



y 



for obtaining the general solution of the system (|T]). Then, we consider 
relations among the elements of A^ to obtain the general solution in the 
form x = A^c of the system (JI} for c ^ 0. This construction haspreviously 
been used by B. Malesevic and B. Radicic [3| (see also |4| and [5|). At the 
end of this paper we will give an application of this results to the matrix 
equation AXB = C. 



2. The main result 



In this section we indicate how technique of an {l}-inverse may be used 
to obtain the necessary and sufficient condition for an existence of a general 
solution of a non-homogeneous linear system. 



Lemma 2.1. The non-homogeneous linear system ([T]) has a solution if and 
only if the last m — a coordinates of the vector d = Qc are zeros, where 
Q g C mxm is regular matrix such that ([2]) holds. 

Proof: The proof follows immediately from Kroneker-Capelli theorem. We 
provide a new proof of the lemma by using the {l}-inverse of the system 
matrix A. The system ([1]) has a solution if and only if c = AAS^c, see 
R. Penrose p). Since A^> is described by the equations (J3J), it follows that 



AA^ = AP 



\ la 


u ' 


_ V 


w _ 



Q = Q 



-i 



\h 


u ' 









Q. 



Hence, we have the following equivalences 



AA^c 



(I-AA^)c=0 



Q- 1 
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-u ' 





In—a 



Qc 



Q~ X Q- 


n- 1 


\ la 


u ' 


" 






-u 




r' - 
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n—a 





Q)c = 
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-u ' 





In— a 



■Ud 



n—a 
I 



d =0 

^n—a w - 



Furthermore, we conclude c 
Theorem 2.2. The vector 



AA^t 



□ 



X 



A^c+{I-A^A)y, 

y g C nxl is an arbitrary column, is the general solution of the system (jTJ), if 
and only if the {1} -inverse A^ of the system matrix A has the form ([3]) for 
arbitrary matrices U and W and the rows of the matrix V{d a — y' a ) + y', n _ a ) 



are free parameters, where Qc 







and P 



y = y 



V'a 



Vn- 



Proof: Since {l}-inverse A^ of the matrix A has the form ([3]), the solution 

of the system x = A^c + (I — A^A)y can be represented in the form 



x 



P 



P 



\ h 


U ' 


_ V 
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w _ 

u ' 


_ V 


w _ 



Qc+i I 
d + ( I 



P 



\ h 


U ' 


_ V 


w _ 



p 



\ L 


u ' 


_ V 


w 



QAjy 
QAPP- 1 



According to Lemma 12.11 and from (j2]) we have 



x 
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\c' a ] 
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w _ 
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Furthermore, we obtain 
P 



x 



Vc> 



P 
P 
P 
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- I I-P 
+ (PP- 1 
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where y' = P y. We now conclude 



x 



P 



Vc'„ 



+ 







■Vy'a + V'n-a 



P 



VK - y' a ) + y' n _ a 



Therefore, since matrix P is regular we deduce that P , , , , 

L V « - Va) + V'n-a 

is the general solution of the system ([1]) if and only if the rows of the matrix 
^( c a ~ Va) + v'n-a are n ~ a f ree parameters. □ 



Corollary 2.3. The vector 



x= (I-A {l) A)i 



y g C nxl is an arbitrary column, is the general solution of the homogeneous 
linear system Ax = 0, A E C- mxn , if and only if the {l}-inverse A^ of the 
system matrix A has the form ([3]) for arbitrary matrices U and W and the 
rows of the matrix —Vy' a + y'( n _ a ) are free parameters, where P~ x y — y' — 



Va 



Vn 



Example 2.4. Consider the homogeneous linear system 

X\ + 2X2 + 3^3 = 
4xi + 5x 2 + 6x 3 = 0. 



(4) 



The system matrix is 
For regular matrices 

Q = 



A 



1 2 3 

4 5 6 



1 

-4 1 



and P 



1 2 1 

3 L 

-I -2 

1 



equality (0|) /ioWs. Rohde's general {1} -inverse A^ l > of the system matrix A 
is of the form 

"1 
A w = P 1 Q 

_ Vu v 12 _ 

According to Corollary \2.3\ the general solution of the system Op is of the 
form 



x = P 









' 











. -Ull 


-V\2 
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p- 1 



2/i 

2/2 
2/3 



where 



P- 1 



Therefore, we obtain 



x 








' 


p 










_ -vu -v 12 


1 


p 
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2 3 




-3 -6 
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2/i + 2y 2 + 3y 3 




-3y 2 - 6y 3 




2/3 













// we take a = —V\\y\ — (2vu — 3vi 2 y 2 — (3t> n — 6v± 2 — l)y 3 as a parameter 
we get the general solution 
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Corollary 2.5. The vector 



x 



A^t 



is the general solution of the system ([I]), if and only if the {1} -inverse A^ 
of the system matrix A has the form ([3]) for arbitrary matrices U and W and 

the rows of the matrix Vc' a are free parameters, where Qc = d = 







Remark 2.6. Similar result can be found in paper B. Males evic and B. 
Radicic [3] . 



Example 2.7. Consider the non-homogeneous linear system 

x 1 + 2x 2 + 3x 3 = 7 
4xi + 5x 2 + 6x 3 = 8. 



(5) 



According to Corollary \2.5\ the general solution of the system |3]) is of the 
form 



x 



P 



1 







' 7 " 
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1 
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= P 
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U12 _ 






7vu - 20-ui2 



// we take a = 7vu — 20f 12 as a parameter we obtain the general solution of 
the system 



x 



P 
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-20 
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We are now concerned with the matrix equation 

AX = C, 
where A e C mxn , X e C nxk and C e C mxk . 



(6) 



Lemma 2.8. The matrix equation ([6]) has a solution if and only if the last 
m — a rows of the matrix C = QC are zeros, where Q G C mxm is regular 
matrix such that (121) holds. 



Proof: If we write X = [X±x Xj,2 • • 
we can observe the matrix equation 



X ik ] and C = [C^ C i2 ■■■ C ik \, then 
3) as the system of matrix equations 



AX ix = C ix 

AX± 2 = C4.2 



AX±k — C^k- 

Each of the matrix equation AX^ = C^, 1 < i < k, by Lemma 12.11 has 
solution if and only if the last m — a coordinates of the vector C'u = QC^ are 
zeros. Thus, the previous system has solution if and only if the last m — a 
rows of the matrix C = QC are zeros, which establishes that the matrix 
equation ([H]) has solution if and only if all entries of the last m — a rows of 
the matrix C are zeros. □ 

Theorem 2.9. The matrix 

X = A^C + (/ - A&A)Y e C nxfc , 

Y e C nxfc is an arbitrary matrix, is the general solution of the matrix equation 
06]) if and only if the {1} -inverse A^ 1 ' of the system matrix A has the form 
03]) for arbitrary matrices U and W and the entries of the matrix 



V(C' a 



Y a ) + Y(n-a) 



are mutually independent free parameters, where QC 
T Y' 

p-ly = Y i 



a 



a 



Yi 



Proof: Applying the Theorem l2.2l on the each system AX^ 
we obtain that 



and 



C iu 1 < i < k, 



X 



i> 



P 



a 



a^i 



I 4' a\i) ' n—a\i . 

is the general solution of the system if and only if the rows of the matrix 
V{C a \.i~ ^aij+^n-aii are n ~ a f ree parameters. Assembling these individual 
solutions together we get that 

a 



X = P 



V(C' a - Yi) + YI 



is the general solution of the matrix equation §§§ if and only if entries of the 



matrix V(C' a -Y>) + Y^ a are 



n — a)k mutually independent free parameters. 

□ 



From now on we proceed with the study of the non-homogeneous linear 
system of the form 

xB = d, (7) 

where B is an n x m matrix over the field C of rank b and d is an 1 x m 
matrix over C. Let R G C nxn and S G C mxm be regular matrices such that 



\ h 


" 









RBS = Eh 



An {l}-inverse of the matrix B can be represented in the Rohde's form 



B™ = S 



\ h 


M ' 


N 


K 



R 



(8) 



(9) 



where M 



a 



yj> 



N 



Vij] and K = 
sponding dimensions b x (n — b), (to 
mutually independent entries. 



[wij] are arbitrary matrices of corre- 
- b) x b and (to — b) x (n — b) with 



Lemma 2.10. The non-homogeneous linear system (JTJ) has a solution if and 
only if the last m — b elements of the row d' = dS are zeros, where S G C mxm 
is regular matrix such that (jS} holds. 

Proof: By transposing the system ((7|) we obtain system B T x T = a? T and by 
transposing the matrix equation (jSJ) we obtain that S T B T R T = Eh. Accord- 
ing to Lemma 12.11 the system B T x T = a? T has solution if and only if the last 
m — b coordinates of the vector S T d T are zeros, i.e. if and only if the last 
m — b elements of the row d' = dS are zeros. □ 

Theorem 2.11. The row 

x = dB {1) +y(I-BB {1) ), 

y g C lxn is an arbitrary row, is the general solution of the system (J7J), if 
and only if the {1} -inverse B^ of the system matrix B has the form §§§ for 
arbitrary matrices N and K and the columns of the matrix (d' b — y' b ) M + y' n _ b 
are free parameters, where dS = d' = [d' b \ 0] and yR~ x — y' — [y' b \ y' n _ b ] ■ 



8 



Proof: The basic idea of the proof is to transpose the system (J7J) and to 
apply the Theorem 12.21 The {l}-inverse of the matrix B T is equal to a 
transpose of the {l}-inverse of the matrix B. Hence, we have 



[B T )^ = (B^f= [S 



\ h 


M ' 


N 


K 



R 



R 



■T 
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N T ~ 


M T 


K T 



s 1 



We can now proceed analogously to the proof of the Theorem 12.21 to obtain 
that 



x 



R 1 



M T (d' b T-y> b T ) + y> n T _ b 



is the general solution of the system B x = d if and only if the rows of 
the matrix M T {d' h T — y' b T ) + y'n-b are n — b free parameters. Therefore, 



x 



[d' b \(d' b -y' b )M + y' n _ b ]R 



is the general solution of the system ([7]) if and only if the columns of the 
matrix (d b — y b )M + y' n _ b are n — b free parameters. □ 



Analogous corollaries hold for the Theorem 12.111 
We now deal with the matrix equation 

XB = D, 
where X e C kxn , B e C nxm and D e C kxm . 



(10) 



Lemma 2.12. The matrix equation ( TTUj) has a solution if and only if the last 
ra — b columns of the matrix D' = DS are zeros, where S G C mxm is regular 
matrix such that ([8]) holds. 



Theorem 2.13. The matrix 

X = DB {1) + Y(I - BB {1) ' 



ifcxn 



Y G C kxn is an arbitrary matrix, is the general solution of the matrix equation 
f fTUj) if and only if the {1} -inverse B^ of the system matrix B has the form 
<Q for arbitrary matrices N and K and the entries of the matrix 



V(D> b -Y>)M + Y{ n „ b) 

are mutually independent free parameters, where DS 
YRT 1 = Y'= [Yl | YU] ■ 



D> = [D' b | 0] and 



3. An application 

In this section we will briefly sketch properties of the general solution of 
the matrix equation 

AXB = C, (11) 

where A e C mxn , X e C nxk , B e C kxl and C G C mxl . If we denote by Y 
matrix product XB, then the matrix equation ( II ip becomes 



AY = C. 



(12) 



According to the Theorem 12.91 the general solution of the system ( 1T2|) can 
be presented as a product of the matrix P and the matrix which has the 
first a = rank(A) rows same as the matrix QC and the elements of the last 
m — a rows are (m — a)n mutually independent free parameters, P and Q 
are regular matrices such that PAQ = E a . Thus, we are now turning on to 
the system of the form 

XB = D. (13) 

By the Theorem 12.131 we conclude that the general solution of the system 
f fT3|) can be presented as a product of the matrix which has the first b = 
rank(B) columns equal to the first b columns of the matrix DS and the rest 
of the columns have mutually independent free parameters as entries, and 
the matrix R, for regular matrices R and S such that RBS = E^. Therefore, 
the general solution of the system (TTTj) is of the form 



X = P 



G a b 


F ' 


H 


L 



R. 



where G a b is a submatrix of the matrix QCS corresponding to the first a 
rows and the first b columns and the entries of the matrices F, H and L are 
nk — ab free parameters. We will illustrate this on the following example. 

Example 3.1. We consider the matrix equation 

AXB = C, 



where A 



1 
-2 



-2 
4 



1 2 1 

,B= 12 1 umlC 

1 2 1 
take Y = XB, we obtain the system 

ay = e 



2 
-4 



. If we 



10 



It is easy to check that the matrix A is of the rank a = 1 and for matrices 



Q 



1 o 

2 1 
the Theorem 



and P 



1 2 
1 



the equality QAP = E a holds. Based on 



2.9\ the equation AY = C can be rewritten in the system form 



AY, 



J.1 



AY, 



4.2 



AY, 



43 



1 

-2 

2 
-4 

1 
-2 



Combining the Theorem \2. S\ with the equality 



c 1 c 2 



(k 







yields 



Y 



41 



1 

2 1 



P 



1 2 

-2 -4 



1 

-2 



12 1 




*i 



42 



P 



Y, 



43 



2w 



v ' 

a 



2vz u + z 22 

v ' 



for an arbitrary matrix Z = 
tion of the system AY = C is 



v - vz 13 + 


^23 


7 


Zl2 #13 


/ 


— ^22 ^23 





Therefore, the general solu- 



Y = P 



1 2 1 

a /3 7 



From now on, we consider the system 

XB = D 
for 

D = P 



1 2 1 




l + 2a 


2 + 2/? l + 2 7 


a (3 7 




a 


/3 7 



11 



There are regular matrices R 



1 


-2 


-1 





1 











1 



1 
— 1 10 and S = 1 such 
-1 1 

that RBS = Ef, holds. Since the rank of the matrix B is b = 1, according to 
the Lemma \2.12\ all entries of the last two columns of the matrix D' = DS 
are zeros, i.e. we have 7 = a, (3 — la. Hence, we get that the matrix D' is 
l + 2a 



of the form D' 



a 







Applying the Theorem \2.13\ we obtain 



X 



l + 2a (l + 2a-t 11 )mn + t 12 (1 + 2a - t u )m 12 + t 



13 



01 



a 



(a-t 2 i)m 1 i + t 



22 



132 

(a - t 12 )m 12 + t 23 



71 



12 



R. 



for an arbitrary matrix T 
system AXB = C is 

X = 



tn £12 £13 

£21 ^22 £23 



Finally, the solution of the 



1 + 2a -Pi -fa Pi Pi 

a - 71 - 72 7i 72 
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